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Abstract 



In these proceedings, we discuss why Functional Renormalization is an essential tool to treat strongly disordered systems. 
More specifically we treat elastic manifolds in a disordered environment. These are goverened by a disorder distribution, which 
after a finite renomalization becomes non-analytic, thus overcoming the predictions of the seemingly exact dimensional reduc- 
tion. We discuss, how a renormalizable field theory can be constructed, even beyond 2-loop order. We then consider an elastic 
manifold imbedded in N dimensions, and give the exact solution for iV — ► oo. This is compared to predictions of the Gaussian 
replica variational ansatz, using replica symmetry breaking. Finally, the effective action at order 1 /N is reported. 



In these proceedings we consider an elastic manifold in a random potential, as prototype for strongly disordered systems. Since 
for all these systems temperature is irrelevant, we will only treat zero temperature. The kind of systems we have in mind 
are domain walls in dirty magnets, contact lines, charge density waves, vortex lattices, to just mention a few. These results 
were obtained in collaboration with Pierre Le Doussal [1-11]. For lack of space we restrict our discussion to the equilibrium. 
Complementary material, especially for the depinning, can be found in the earlier review [12]. 



The simplest experimental realization is an Ising magnet. Imposing boundary conditions with all spins up at the upper and all 
spins down at the lower boundary (see figure 1), at low temperatures, a domain wall separates a region with spin up from a region 
with spin down. In a pure system at temperature T = 0, this domain wall is completely flat. Disorder can deform the domain 
wall, making it eventually rough again. Figure 1 shows, how the domain wall is described by a displacement field u(x). Another 
example is the contact line of water (or liquid Helium), wetting a rough substrate. A realization with a 2-parameter displacement 
field u(x) is the deformation of a vortex lattice: the position of each vortex is deformed from x to x + u(x). A 3-dimensional 
example are charge density waves. 

All these models have in common, that they are described by a displacement field x 6 R d — > u(x) £ R . For simplicity, 
we set N = 1, if not explicitly stated otherwise. After some initial coarse-graining, the energy Ti = Tt c \ + T~Loo consists out of 
two parts: the elastic energy Ti c \ and the disorder energy Tioo 



1 Introduction 



2 Physical realizations, model and observables 




(2.1) 



We choose the disorder at the microscopic scale Gaussian, with correlations 



V(u, x)V(u',x') := 5 d {x - x')R(u - u) . 



(2.2) 



The most interesting observable is the roughness-exponent from the behavior of the correlation function 



[u(x) - u{y)f ~ \x-y\ 



(2.3) 



Other observables are higher correlation functions or the free energy. 
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Figure 1: An Ising magnet at low temperatures forms a domain wall described by a function u(x) (right). An experiment on a thin Cobalt film (left) [13]; with 
kind permission of the authors. 

3 Dimensional reduction 

There is a beautiful and rather mind-boggling theorem relating disordered systems to pure systems (i.e. without disorder), which 
applies to a large class of systems, e.g. random field systems and elastic manifolds in disorder. It is called dimensional reduction 
and reads as follows [14]: 

Theorem: A d-dimensional disordered system at zero temperature is equivalent to all orders in perturbation theory to a pure 
system in d — 2 dimensions at finite temperature. 

Let me give an example: The thermal expectation value for the 2-point function scales as {[u{x) — u(y)] 2 ) ~ |x| 2 ~ d . 
Making the dimensional shift implied by dimensional reduction implies that the disorder-averaged 2-point function at zero 
temperature is 

4— d 



[u(x)-u(0)Y 



C = 



d 



(3.1) 



We will see later that this is not true; but remains an important benchmark due to fact that the "theorem" is correct to all orders 
in the disorder strength and its moments (i.e. when expanding in _R"(0), R""(0), a.s.o.). 

4 The Larkin-length 

To understand the failure of dimensional reduction, let us turn to an interesting argument given by Larkin [15]. He considers a 
piece of an elastic manifold of size L. If the disorder has correlation length r, and characteristic potential energy e, this piece 
will typically see a potential energy of strength -Edo = e(-^) d//2 '. On the other hand, there is an elastic energy, which scales like 
-Ed = cL d ~ 2 . These energies are balanced at the Larkin-length L = L c with L c = (%-r d ) 1//< - 4_d '. More important than this 
value is the observation that in all physically interesting dimensions d < 4, and at scales L > L c , the membrane is pinned by 
disorder; whereas on small scales elastic energy dominates. This means that d — 4 is the upper critical dimension. 

5 The functional renormalization group (FRG) 

Let us now discuss a way out of the dilemma, posed by dimensional reduction: We would like to make an e = 4 — d expansion. 
On the other hand, dimensional reduction tells us that the roughness is £ — ^2 ( se e (3.1)). Even though this is systematically 
wrong below four dimensions, it tells us correctly that at the critical dimension d — 4, where disorder is marginally relevant, 
the field u is dimensionless. This means that having identified any relevant or marginal perturbation (as the disorder), we find 



-R"(u) 




-R"(u) A 



renormalization 




Figure 2: Change of —R"(u) under renormalization and formation of the cusp. 
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Figure 3: Generation of the cusp, as explained in the main text. 



immediately another such perturbation by adding more powers of the field. We can thus not restrict ourselves to keeping solely 
the first moments of the disorder, but have to keep the whole disorder-distribution function R(u). Thus we need afunctional 
renormalization group treatment (FRG). Functional renormalization is an old idea going back to the seventies, and can e.g. be 
found in [16], by Wegner and Houghton. For disordered systems, it was first proposed in 1986 by D. Fisher [17]. Performing an 
infinitesimal renormalization, i.e. integrating over a momentum shell a la Wilson, leads to the flow diR(u), with (e = 4 — d) 

d e R{u) = (e - 4C) R(u) + (uR'(u) + \R" (u) 2 - R"(u)R"(0) . (5.1) 

The first two terms come from the rescaling of R and u respectively. The last two terms are the result of the 1-loop calculations, 
see e.g. [4] . 

More important than the form of this equation is its actual solution, sketched in figure 2. After some finite renormalization, 
the second derivative of the disorder R"(u) acquires a cusp at u = 0; the length at which this happens is the Larkin-length. How 
does this overcome dimensional reduction? To understand this, it is interesting to study the flow of the second and forth moment. 
Taking derivatives of (5.1) w.r.t. u and setting u to 0, we obtain 

deR"(0) = (e-2()R"(0) + R"'(0) 2 —> {e-2QR"(0) (5.2) 
d e R""(0) = £ J R""(0) + 3 J R""(0) 2 + 4 J R'"(0) J R'""(0) — > eR"" (0) + 3R"" (0) 2 . (5.3) 

Since R(u) is an even function, R"'(0) and R""'(0) are 0, which we have already indicated in Eqs. (5.2) and (5.3) . The 
above equations for R"(0) and R""(0) are in fact closed. Equation (5.2) tells us that the flow of R"(0) is trivial and that 
£ = e/2 = — This is exactly the result predicted by dimensional reduction. The appearance of the cusp can be inferred 

from equation (5.3). Its solution is R (0)1/ = j 2 -, — r— , with c := _R""(0)U=o. Thus after a finite renormalization 

1 — 3 cl e e * — 1 j/e 

R""(0) becomes infinite: The cusp appears. By analyzing the solution of the flow-equation (5.1), one also finds that beyond 
the Larkin-length R"(0) is no longer given by (5.2) with R"'{0) 2 = 0, but R"'(0) 2 -> R"'{0 + ) 2 = lim u ^o R"'{u) 2 , which 
is non-zero after the cusp. Renormalization of the whole function thus overcomes dimensional reduction. The appearance of 
the cusp also explains why dimensional reduction breaks down. The simplest way to see this is by redoing the proof for elastic 
manifolds in disorder, which in the absence of disorder is a simple Gaussian theory. Terms contributing to the 2-point function 
involve _R"(0), TR""(0) and higher derivatives of R(u) at u — 0, which all come with higher powers of T. To obtain the limit 
of T — > 0, one sets T = 0, and only R"(0) remains. This is the dimensional reduction result. However we just saw that R""(0) 
becomes infinite. Thus R""(0)T may also contribute, and the proof fails. 

6 The cusp and shocks 

Let us give a simple argument of why a cusp is a physical necessity, and not an artifact. The argument is quite old and appeared 
probably first in the treatment of correlation-functions by shocks in Burgers turbulence. It became popular in [18]. Suppose, we 
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Figure 4: Results for £ in the random bond case. 



want to integrate out a single degree of freedom, whose average position due to the elastic energy connecting it to its neighbors is 
u. This harmonic potential and the disorder term are represented by the parabola and the lowest curve on figure 3(a) respectively; 
their sum is the remaining curve. For a given disorder realization, the minimum of the potential as a function of u is reported on 
figure 3(b). Note that it has non-analytic points, which mark the transition from one minimum to another. Taking the derivative 
of the potential leads to the force in figure 3(c). It is characterized by almost linear pieces, and shocks (i.e. jumps). Calculating 
the force-force correlator, the dominant contribution in its decay for small distances is due to the presence of shocks. Their 
contribution is proportional to their probability, itself proportional to the distance between the two observable points. This leads 
to F(u)F(0) — F(0) 2 — c\u\, with some numerical coefficient c. 

7 Beyond 1 loop? 

Functional renormalization has successfully been applied to a bunch of problems at 1-loop order. From a field theory, we 
however demand more. Namely that it allows for systematic corrections beyond 1-loop order; be renormalizable; and thus 
allows to make universal predictions. However, this has been a puzzle since 1986, and it has even been suggested that the theory 
is not renormalizable due to the appearance of terms of order e 3//2 [19]. Why is the next order so complicated? The reason is that 
it involves terms proportional to R"'(0). A look at figure 3 explains the puzzle. Shall we use the symmetry of R(u) to conclude 
that R'" (0) is 0? Or shall we take the left-hand or right-hand derivatives, related by 

R"'(0 + ) := lim R"'(u) = - lim R"'{u) =: -R"'(Q-). (7.1) 

u^O u — 

In the following, I will present the solution of this puzzle, at 2-loop order and large N. The latter approach allows for another 
independent control-parameter, and sheds further light on the cusp-formation. 

8 Results at 2-loop order 

For the flow-equation at 2-loop order, the result is [1,4,20,21] 

diR(u) = (e~4()R(u) + (uR'(u) + ±R"(u) 2 -R"(u)R"(0) 

±i (R"(u) - R"(0)) R'"{u) 2 - \R'"(Q+fR"[u) . (8.1) 

The first line is the result at 1-loop order, already given in (5.1). The second line is new. The most interesting term is the 
last one, which involves _R"'(0 + ) and which we therefore call anomalous. The hard task is to fix the prefactor (— i). We 
have found five different prescriptions to calculate it: The sloop-algorithm, recursive construction, reparametrization invariance, 
renormalizability, and potentiality [1,22]. For lack of space, we restrain our discussion to the last two ones. At 2-loop order the 
following diagram appears 

A R" 

,-AA-- — ► - (R"(u) - R"(0)) R"'(uf - -R"(u)R"'(0 + ) 2 (8.2) 
«""!< 2 2 



leading to the anomalous term. The integral (not written here) contains a subdivergence, which is indicated by the box. Renor- 
malizability demands that its leading divergence (which is of order 1/e 2 ) be canceled by a 1-loop counter-term. The latter is 
unique thus fixing the prefactor of the anomalous term. 

Another very physical demand is that the problem remain potential, i.e. that forces still derive from a potential. The force- 
force correlation function being —R"(u), this means that the flow of R'(0) has to be strictly 0. From (8) one can check that this 
does not remain true if one changes the prefactor of the last term in (8); thus fixing it. 

Let us give some results for random-bond disorder (short-ranged potential-potential correlation function). For this, we have 
to solve (8.1) numerically, with the result £ = 0.20829804e ± 0.006858e 2 . This compares well with numerical simulations, see 
figure 4. 
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Figure 5: Results for the roughness £ at 1- and 2-loop order, as a function of the number of components TV. 



9 Finite N 

Up to now, we have studied the functional RG in two cases: For one component TV = 1 and in the limit of a large number of 
components, TV — > oo. The general case of finite TV is more difficult to handle, since derivatives of the renormalized disorder 
now depend on the direction, in which this derivative are taken. Define amplitude u := \u\ and direction u := u/\u\ of the 
field. Then deriving the latter variable leads to terms proportional to 1/u, which are diverging in the limit of u — ► 0. This poses 
additional problems in the calculation, beyond the case TV = 1. At 1-loop order everything is well-defined [19]. We have found 
a consistent RG-equation at 2-loop order (see [12] and unpublished): 

The fixed point equation has to be integrated numerically, order by order in e. The result, specialized to directed polymers, 
i.e. e = 3 is plotted on figure 5. We see that the 2-loop corrections are rather big at large TV, so some doubt on the applicability 
of the latter down to e = 3 is advised. However both 1- and 2-loop results reproduce well the two known points on the curve: 
£ = 2/3 for TV = 1 and £ = for TV = oo. The latter result has been given in section 10 Via the equivalence [24] of the 
directed polymer problem in TV dimensions treated here and the KPZ-equation of non-linear surface growth in TV dimensions, 
we conclude that d ~ 2.4 is the upper critical dimension of KPZ. 

10 Large N 

In the last section, we have discussed renormalization in a loop expansion, i.e. expansion in e. In order to independently check 
consistency it is good to have a non-perturbative approach. This is achieved by the large-TV limit, which can be solved analytically 
and to which we turn now. We start from the disorder-averaged energy with disorder correlator B(u 2 ) = R(\u\) where we use 
an TV-component field u. We then calculate the free energy in presence of a source j, and finally the effective action T(u) via a 
Legendre transform. For large TV the saddle point equation reads [2] 

B'(u 2 ab ) = 73' (Xab) , Xab = u ab + 2Th + Ah[B'{u 2 ab ) ~ S'(0)] (10.1) 

This equation gives the derivative of the effective (renormalized) disorder 73 as a function of the (constant) background field 
u ab — i u a — Ub) 2 in terms of: the derivative of the microscopic (bare) disorder 73, the temperature T and the integrals I n := 

Ik (fe2 +m 2)™ • 

The saddle-point equation can be turned into a closed functional renormalization group equation for B by taking the derivative 
w.r.t. m (restricting ourselves to T = 0): 

8,B(x) = ~^B(x) = (e - 4QB[x) + 2(xB'(x) + \& {xf - B'(x)B'(0) (10.2) 

This is a complicated nonlinear partial differential equation. It is therefore surprising, that one can find an analytic solution. 
(The trick is to write down the flow-equation for the inverse function of B'(x), which is linear.) Let us only give the results 
of this analytic solution: First of all, for long-range correlated disorder of the form B'(x) ~ x^ 1 , the exponent £ can be 
calculated analytically as f = 2 (i+ 7 ) • ^ a 8 rees wim the replica-treatment in [25] and the 1-loop treatment in [19]. For short- 
range correlated disorder, Q = 0. Second, it demonstrates that before the Larkin-length, B(x) is analytic and thus dimensional 
reduction holds. Beyond the Larkin length, 73" (0) = oo, a cusp appears and dimensional reduction is incorrect. This shows 
again that the cusp is not an artifact of the perturbative expansion, but an important property even of the exact solution of the 
problem (here in the limit of large TV). 
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11 Relation to Replica Symmetry Breaking (RSB) 



There is another treatment of the limit of large N given by Mezard and Parisi [25]. They make a Gaussian variational ansatz of 
the form 



1 - f 1 

ftg[w] = I Ua(x)(-V 2 +m 2 )u a (x)- — ^ a ab u a (x)u b (x) 



2T 2 



(11.1) 



which becomes exact for N — » oo. The art is to make an appropriate ansatz for a ab . The simplest possibility, a a t = a for 
all a 7^ b reproduces the dimensional reduction result, which breaks down at the Larkin length. Beyond that scale, a replica 



symmetry broken (RSB) ansatz for a a t is necessary, of the form o~ a b = 




Parisi has shown that this infinitely often 



replica- symmetry broken matrix can be parameterized by a function [o"](z) with z G [0, 1] where z = describes distant states, 
whereas z = 1 describes nearby states. The solution of the large- N saddle-point equations leads to the curve depicted in figure 
11 Knowing it, the 2-point function is given by (u k u- k ) = pnp^r U + Jo 1 ■ 

What is the relation between the two approaches, which both pretend to calculate the same 2-point function? Comparing the 
analytical solutions, we find that the 2-point function given by FRG is the same as that of RSB, if in the latter expression we only 
take into account the contribution from the most distant states, i.e. those for z between and z m (see figure 6). To understand 
why this is so, we have to remember that the two calculations were done under quite different assumptions: In contrast to the 
RSB-calculation, the FRG-approach calculated the partition function in presence of an external field j, which was then used to 
give via a Legendre transformation the effective action. Even if the field j is finally turned to 0, the system will remember its 
preparation, as is the case for a magnet. 

By explicitly breaking the replica-symmetry through an applied field, all replicas will settle in distant states, and the close 
states from the Parisi-function [a] (z) + m 2 (which describes spontaneous RSB) will not contribute. However, we found that the 
full RSB-result can be reconstructed by remarking that the part of the curve between z m and z c is independent of the infrared 
cutoff 77i, and then integrating over m [2] (m c is the mass corresponding to z c ): 



[u k u^ k 



+ 



dR?„(0) 



1 



1 



(11.2) 



We also note that a similar effective action has been proposed in [18]. 
correct FRG 2-point function, as it should. 



While it agrees qualitatively, it does not reproduce the 



In a graphical notation, we find [11] 

6BW = 



12 Corrections at order 1/N 




O = B"(xat)(l-4A d I 2 (p)B"( Xab )) 



B(Xab) 



(12.1) 
(12.2) 



[a](z)+m 



from UV-cutoff 



IR-cutoff 




Zm 5= 1 

Figure 6: The function [a] (it) + m 2 as given in [25]. 



6 



where the explicit expressions are given in [11]. By varying the IR-regulator, one can derive a /3-function at order 1/N, see 
[11]. At T = 0, it is UV-convegent, and should allow to find a fixed point. We have been able to do this at order e, showing 
consistency with the 1-loop result, see section 9. Other dimensions are more complicated. 

A /3-function can also be defined at finite T. However since temperature is an irrelevant variable, it makes the theory non- 
renormalizable, i.e. in otder to define it, one must keep an explicit infrared cutoff. These problems will be treated in a forthcoming 
publication. 

13 Perspectives 

Other interesting problems have been treated by the above methods, especially dynamic problems (see [12] for a review); and 
many more are now in reach. Some open points have already been raised in these notes, others are the strong disorder phase of 
random field problems, or wether FRG can also be applied to spin-glasses. We have to leave these problems for future research 
and as a challenge for the reader to plunge deeper into the mysteries of functional renormalization. 
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